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ABSTRACT 

We have discovered an intermediate polar (IP) in the 100 ks Chandra observation of Baade's Window (BW), 
a low extinction region at about 4° south of the Galactic Center. The source exhibits large X-ray modulations at 
a period of 1028.4 s in the 0.3-8 keV band. The X-ray spectral fit with a power law model shows the integrated 
spectrum is intrinsically hard (photon index, Y = 0.44 ± 0.05) and moderately absorbed (A/h = 1 .5 ± 1 .0 x 10 21 
cm" 2 ). The relatively poor statistics only allows for a mild constraint on the presence of an iron emission line 
(equivalent width = 0.5 ±0.3 keV at 6.7 keV). Quantile analysis reveals that the modulations in the X-ray flux 
strongly correlate with spectral changes that are dominated by varying internal absorption. The X-ray spectrum 
of the source is heavily absorbed (Nh > 10 22 cm -2 ) during the faint phases, while the absorption is consistent 
with the field value (~ 10 21 cm" 2 ) during the bright phases. These X-ray properties are typical signatures of 
IPs. Images taken with the IMACS camera on the Magellan 6.5m telescope show a faint (V ~ 22), relatively 
blue object (Bq — Vq > 0.05) within the 2a error circle of the Chandra source, which is a good candidate for 
being the optical counterpart. If we assume a nominal range of absolute V magnitude for a cataclysmic variable 
(My ~ 5.5- 10.5) and the known reddening in the region (Ay = 1.4 at > 3 kpc), the source would likely be at 
a distance of 2-10 kpc and not in the local solar neighborhood. The corresponding average X-ray luminosity 
would be 6 x 10 31 - 10 33 ergs s" 1 in the 2-8 keV band. Assuming the space density of IPs follows the stellar 
distribution, which is highly concentrated in the Galactic Bulge, the source is probably a relatively bright IP 
(~ 10 33 ergs s" 1 if it is at 8 kpc) belonging to the Galactic Bulge X-ray population, the majority of which is 
now believed to be magnetic cataclysmic variables. 

Subject headings: Galaxy: bulge — X-ray: binaries — cataclysmic variables 



1. INTRODUCTION 

Intermediate Polars (IPs) are a type of magnetic cata- 
clysmic variable (CV), where the magnetic field of the ac- 
cretor, a white dwarf (WD), disrupts the inner portion of the 
accretion disk and ch annels the accre tion flow into the mag- 
netic poles of the WD (lPattersonll998l) . As the WD spins, this 
channeling gives rise to a pulsation, a tell sign of an IP, which 
is typically found in the period range of ~ 100 - 1000s. The 
magnetic field on the surface of the WD in IPs is typically 
about 10 6 ~ 7 G, and the ratio of the spin to orbital period of 
IPs is found in a wide range from 0.01 to <1. Some IPs are 
on the evolutionary path to their cousins, the polars, where 
the orbital and WD spin periods are usually locked under the 
strong magne tic fields that convert the whole acc retion flow 
into a stream (iNorton. Wynn & Somersc ales 2004). Magnetic 
CVs provide unique astrophysical laboratories for studying 
physics in extreme conditions, stellar and binary evolution, 
which is important to make a census of Galactic population 
and understand their evolution. There are about 33 confirmed 
and 72 candidate IPs today, and the full catalogue is found in 
Ritte r&Korbl(l2003l) 3 . 

The 72 IP candidates contain seven periodic X-ray sources 
at the Galactic Center (GC), which are part of ~ 2300 X-ray 
sources disc overed in the deep Chandra observations of the 
Sgr A* field (Muno et al. 2003a). The exact nature of the ma- 
jority of the GC X-ray population still remains elusive: direct 
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identification at other wavelength is difficult due to high ob- 
scuration by dust and source confusion due to high star den- 
sity. Therefore, the discovery of these IP candidates is partic- 
ularly interesting since it supports the idea of magnetic CVs 
as the le ading candidates for the majority of t he GC X-ray 
sources dMuno et al]l2004t ILaycock et al.ll2005l) . If true, this 
also implies the GC X-ray population is very old. 

In order to explore the X-ray population in the Galactic 
Bulge (GB) without obscuration by the intervening dust, we 
have observed three low extinction Windo ws within 4° of the 
GC with the Chandra X-ray observatory ( Ivan den Berg et ail 
l2006ll2009l:rHong et al.l2009b . The initial population study on 
the X-ray sources in these Windows and four GB fields has re- 
vealed that the GB X-ray sources extend out to > 1 .4° from 
the GC with a projected source density that follo ws a 1 /0 re- 
lation , where 8 is the angular offset from the GC dHong et all 
2009). In addition, the similarity of the X-ray spectra of the 
hard X-ray sources in these fields - an intrinsically hard con- 
tinuum with the presence of an iron emission line - indicates 
a single class of sources, likely magnetic CVs, could make up 
both the GC and GB X-ray populations. 

As part of our efforts to identify the nature of these sources, 
we have searched for periodic modulations in the X-ray emis- 
sion of the bright X-ray sources (net counts > 250 in the 0.3-8 
keV ran ge) in the Win dow fields using the Lomb-Scargle al- 
gorithm (Scargle 1983). As a result, we have discovered an 
IP in Baade's Window (BW) based on a strong pulsation cor- 
related with the X-ray spectral variation (©. We explore the 
properties of this X-ray source (® and the potential optical 
counterpart (® that lead to its identification as an IP. We dis- 
cuss the possibility that the system belongs to the GB X-ray 
"opulation (©. 
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2. OBSERVATION AND TIMING ANALYSIS 
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FIG. 1.— The raw X-ray sky image around CXOPS J180354.3-300005. 
The (red) circle indicates the source region used for the aperture photometry. 
The morphology of the event distribution is consistent with the expected point 
spread function at the source location. 

We have observed BW on 2003 July 9 (Obs. ID 3780), 
Stanek's Window (SW) on 2004 February 14/15 (Obs. ID 
4547 and 5303), and Limiting Window (LW) on 2005 August 
19/22 and October 25 (Obs. I D 5934, 6362 and 6365) with the 
Chandra ACIS-I instrument (van den Berg et afll2006l 120091 
iHong et al.l2009l) . The data were analyzed as a part of our sur- 
vey program, the Chandra Multi-wavelength Plane (ChaM- 
Plane) survey, which is designed to constrain the Galactic 
populat ion of low luminosity accretion sources, CVs in par- 
ticular (Grindl ay et al.l 120051) . In summary, we search for 
X-ray point sources using a wavelet algorithm (wavdetect, 
Freem an" et alj |2~002) and perform aperture photometry to ex- 
tract the basic X-ray prope rties. The details of the analysis 
procedures are described in Hong et al. (2005, 2009]) . 

In order to find periodic X-ray modulations in the light 
curves of these X-ray sources, we have performed a Lomb- 
Scargle periodicity search on the bright X-ray sources (net 
counts > 250 in 0.3-8.0 keV) discovered in the 100 ks obser- 
vation of the Window fields - BW (Obs. ID 3780), Stanek 
Window (Obs. ID 4547 and 5303) and Limiti ng Window 
(Obs. ID 5934, 6362 and 6365) (IHong et alJ2009l) . The events 
are selected by the good time intervals (GTIs), where the 
background fluc tuation is less than 3cr above the mean level 
(Hong et al. 2005). The arrival time of each photon is bary- 
center corrected using the CI AO tool axbary 4 . For each 
source, we generate the light curves in time bins of four multi- 
ples (1, 4, 8 and 16) of the CCD integration time (3.2 s). Then 
we apply the Lomb-Scargle algorithm and calculate the power 
spectrum for the light curves at all four time resolutions. 

We have found that one source in BW, CXOPS J 180354.3- 
300005 exhibits a clear sign of a periodic modulation with 
greater than 99% of the confidence level (CL). The source was 
detected with 5 10 ± 24 net counts in the broad band (Bx, 0.3- 
8 keV) after background subtraction. Fig. 1 shows the raw 
sky image around the source marked by the (red) circle of the 
95% point spread function (PSF). The source only appears to 
be extended due to the large offset (8.5') from the aimpoint of 
the instrument, but a simulated PSF by a CIAO tool mkpsf 
at the source position closely resembles the event distribution 
of the source, indicating that the source is consistent with a 
point source. 

Fig. 2a shows the power spectrum of the source in the fre- 



quency domain along with horizontal lines indicating 90, 95 
and 99% CLs. This power spectrum is based on the light curve 
at the 12.8 s time resolution and the results based on the other 
three time resolutions are very similar. The power spectrum 
indicates the source exhibits pulsations at a period of 1028.4 
± 3.8 s. The error of the period estimate is the lcr equivalent 
spread of the peak in the power spectrum, which is calculated 
by a fit to the peak with a Gaussian function in the period 
domain. The observed period is somewhat close to a possi- 
ble spurious period (1000 s). Sources falling near the edges 
or node boundaries of the chips can exhibit spurious modu- 
lations at a period of 707 or 1000 s (or their harmonics) due 
to the dither motion of the Chandra X-ray observatory. How- 
ever, we consider the periodic modulation found in CXOPS 
J 180354.3-300005 to be real because first, the source is not 
near the chip edges or the node boundaries and second, the 
X-ray spectral properties of the source strongly correlate with 
the flux modulation as demonstrated below. 

Fig. 2b shows the folded light curves at the pulsation period 
in the soft (0.3-2.0 keV), hard (2.0-8.0 keV) and broad (0.3- 
8.0 keV) bands. The folded light curve of the broad band is 
shifted up by 1 cts ks" 1 for clarity. The smooth fines are cal- 
culated using the LOWESS algorithm (IClevelandlll994h and 
color-coded by the phase for a later reference in the spectral 
analysis (§3). The data points are calculated at the 10 equal- 
sized phase bins of width 0. 1 . The net count per bin ranges 
from ~ 18 to 90 in the 0.3-8.0 keV range. 

We define the modulation depth as (R mll x-Rmm)/(Rmax + 
^min) where /? max and R m i n are the maximum and minimum of 
the fitted amplitudes respectively. Using a fit to the 10 data 
points with a sinusoidal function, we estimate the modulation 
depth is (90 ± 10)%, (42 ± 10)% and (53 ± 8)% for the soft, 
hard, and broad bands respectively. The modulation ampli- 
tude varies with the energy band, which we will explore in 
more detail using quantile analysis in §|3] 

B W was also observed with the EPIC cameras on the XMM- 
Newton observatory in two separate pointings with ~ 20-25 
ks exposure each on 2002 March 1 1 and 2004 September 30. 
The data from both observations are publicly available. Our 
source was detected as 2XMM J180354.3-300004 5 at large 
offsets (5.7' and 7.6 ') from the aimpoint in both observations 
d Watson et a l. 2009). Unfortunately, in the four datasets out 
of six total (three cameras and two pointings), the source fell 
near a chip gap. In both PN observations, which would have 
provided an interesting result because of the superior sensitiv- 
ity at high energies than the Chandra ACIS instruments, the 
central part of the PSF overlapped with a chip gap, rendering 
the data practically unusable. Three MOS observations (in- 
cluding one with the source marginally close (~30") to the 
chip edge) appear to produce relatively clean dataset of the 
source, but the poor statistics (net counts: ^60-80 each in 
0.3-8 keV) did not allow for any significant detection of pul- 
sation. Therefore, in the following, we mainly consider the 
Chandra data for the X-ray observation of the source. See 
Table 1 for the overall flux estimates of the source from the 
XMM-Newton observations, which is consistent with the same 
from the Chandra data. 

3. X-RAY SPECTRUM AND VARIATION 

Table 1 summarizes the results of the spectral model fits to 
the integrated X-ray spectrum of the source using a power law 
model with and without an iron emission line at 6.7 keV. The 
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FIG. 2.— The power spectrum (a) of CXOPS J18O354.3-30O0O5 in the broad band (0.3-8.0 keV) and the folded light curves (b) in the soft (0.3-2.0 keV), hard 
(2.0-8.0 keV), and broad bands (phase = at JD 2452830.0). The primary modulation is found at 9.727 ± 0.036 X 10 -4 Hz or 1028.4 ± 3.8 s. The horizontal 
lines in (a) represent the confidence levels (CLs) of 90, 95 and 99%. The smooth lines in (b) are generated by the LOWESS alghorithm and color-coded by the 
phase for a later reference (Fig. [4). The data points are calculated in the equal phase bins (0.1). The modulation depth is estimated to be (90± 10)%, (42± 10)% 
and (53 ± 8)% for the soft, hard, and broad bands respectively. The folded light curves in the broad band are shifted up by 1 cts ks~' for clarity. 



TABLE 1 

Spectral Model Fits and Quantile Analysis of the integrated X-ray spectrum of CXOPS J 180354. 3-300005 



Spectral Fits Quantile Analysis Flux 

Model F EW p7DoF~ %22 F ™~ 0.5 - 2 keV 2 - 8 keV 

(xl0 22 cnr 2 ) (keV) (x!0 22 cnT 2 ) (keV) (xlO~ 13 ergs cm' 2 s" 1 ) 

Chandra 

PL 0.15(10) 0.44(5) - 12.6/22 0.28(21) 0.49(18) - 0.13(1) 1.20(7) 

PL+Fe 0.22(10) 0.53(5) 0.5(3) 10.3/20 0.39(22) 0.63(18) 05 0.14(1) 1.19(7) 

XMM-Newton MOS 

PL 0.6l!gj 4 0-60$39 : 0.19(12) 1.25(43) 

For the Chandra data, the flux estimates are based on the spectral model fits. The errors are of statistical origin, based on the photon counts using small-number 
statistics by Gehrels et al. 1 1986). The PL+Fe model in quantile analysis assumes a 0.5 keV EW for the iron line at 6.7 keV. In the XMM-Newton data, we 
combine three datasets and use quantile analysys under a simple power law. The datasets comprise the MOS 1 observation in the 2002 pointing, and both the 
MOS 1 & MOS2 observations in 2004. The flux estimate is based on the quantile analysis of the combined dataset, and the errors are the quadratic sum of the 
statistical error and the standard deviation of the t hree datasets. The r eported flux in 0.2-12 keV is 1 .9(2) X 10~ 13 and 1 .6(2) x 10~ 13 ergs cm~ 2 s -1 for the 2002 
and 2004 XMM-Newton observations respectively (Watson et al. 2009, see also http : / / x mm .esac.esa. int /xsa/). DoF stands for degrees of freedom. 



Fe He-a XXV line is considered because many of the GC X- 
ray sources exhibit this li ne more prominently than the neu- 
tral Fe line at 6.4 keV dWang et al.ll2002t iMuno et al.l 120031; 
iHong et ai]|2009l) . Both models produce acceptable fits, and 
the spectrum is intrinsically hard (T ~ 0.4) and moderately 
absorbed (Nmi ~ 0.2). Fig. [3] shows the spectral model fit 
to the integrated spectrum of the source using the power law 
plus the iron line. There is a hint of an emission line in the 
6-7 keV energy range in the spectrum but the poor statis- 
tics does not provide a significant constraint on the presence 
of the iron line. The estimated equivalent width (EW) of the 
iron line at 6.7 keV is 0.5 ± 0.3 keV. In the case of the ther- 
mal Bremsstrahlung or thermal plasma model, we can only 
assign a lower limit for the temperature (> 10 keV with 95% 
confidence level for thermal Bremsstrahlung or > 7.0 keV for 
MeKaL) and we cannot confirm this apparent high tempera- 
ture is intrinsic due to low statistics and lack of the spectral 
information at high energies (> 10 keV). 

The total column density in the direction of CXOPS 
J 180354.3-300005 is estimated t o be N H2 2 = 0.25(8) by 
iMarshall et all (120061) or ^25(^)^ [Surnil (120041) . using N H2 2 = 
0.179 Ay dPredehl & Schmitll 1995b . The former estimates are 
given in a set of distances starting from 2.75 kpc for the BW 
field and the latter estimates are given at a finer angular reso- 
lution (~ 0.5' - 2.0') than the former (~ 8.0'). There are also 
ot her estimates for the e xtinction in t he BW such as ~ 0.3 4 
bv lSchlegel etail d!998l) or ~ 0.41 bv lDrimmel eTail d2003l) . 
but the underlying models of these estimates are not accurate 
in this region. 



In order to explore the spectral variation as a function of the 
modulation phase, we use quantile analysis. Quantile analysis 
is a bias-free spectral classification method, suitable for faint 
sources or for investigating subsets of the data th at do not al- 
low f or spectral model fits due to low statistics dHong et al.l 
2004). Unlike the conventional hardness ratio, which is sub- 
ject to the spectral bias intrinsic to the choice of the sub- 
energy ranges, quantile analysis provides a reliable spectral 
measure of sources even with counts as low as ~ 10 because 
it takes full advantange of th e given statistics without sub di- 
viding the energy range (e.g. lvan den Berg et al.l l2~006). 

We calculate the quantile values of the X-ray spectra of 
the source as a function of modulation phase using a slid- 
ing phase window of a fixed size (0.1). Fig. [4] shows such 
a quantile diagram. The energy quantile E p corresponds to 
the energy below which p% of the counts are detected in 
the energy range (0.3-8.0 keV in Fig. 0). For instance, £50 
means the median energy, and £25 and £75 are two quartiles. 
Note the definitio n of the x-ax i s in F ig.[4]is different from the 
one suggested by Hong et al. (2004). The new definition is 
log 10 (£ 50 /£io)/log 10 (£hi/£'io) in general, where £i and £ hi 
are the lower and upper bound of the energy range respec- 
tively (£i o =0.3 and £ hi =8.0 keV in Fig. g]). We believe the 
new definition is more reflective of the instrument response 
and the confined range (0-1) allows for statistically more uni- 
form response throu ghout the phase space of quantile diagram 
dHong et alj|2009bl) . 

For easy interpretation, the quantile diagram is overlaid 
with power law and thermal Bremsstrahlung model grids. The 
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FIG. 3. — The spectral model fit to the integrated X-ray spectrum of CXOPS 
J180354. 3-300005 under a power law plus an iron line at 6.7 keV. See also 
Table 1 . The spectral fit shows a hint of the iron emission line. 

filled circles are color-coded to match the phase of the smooth 
line in the folded light curve in Fig. 3b for easy comparison 
and the data points with the error bars are from the same 10 
data points in the folded light curve. Table 1 also shows the 
spectral parameters estimated by quantile analysis of the inte- 
grated spectrum (marked by a hollow cross in Fig. |4j and the 
results are consistent with the spectral fit. 

The quantile diagram reveals dramatic spectral changes as a 
function of pulsation phase, illustrated by a large loop formed 
counter-clockwise as the X-ray flux modulates in a full cy- 
cle (see the track in Fig. ®>). The spectral variation is ex- 
pected from the energy dependence of the modulation depth 
seen in Fig. [2] In the quantile diagram, we can also see that the 
spectral changes are dominated by varying extinction. During 
the bright phases, the spectrum is relatively unabsorbed with 
Nhzi % 1) l ess than or consistent with the field extinction (e.g. 
the yellow and red points in Fig. I2J5 & ^p), but during the 
faint phases, especially around phases 0.3 - 0.6 (e.g. the light 
blue points), the spectrum is heavily absorbed with Nmi > 5. 
Since the field extinction is very low, the additional increase 
in extinction must be intrinsic to the source. Note the long 
axis of the data loop in the quantile diagram is largely parallel 
to the direction related to changes in absorption for both spec- 
tral models as shown in Fig. |4] In fact, this is true in general 
regardless of the specral models. 

Table Q] also shows the results of the XMM-Newton data 
for comparison, for which we combine three relatively clean 
datasets (the MOS1 observation in 2002, both MOS1 and 
MOS2 in 2004). The overall spectral properties and the flux 
estimates of the XMM-Newton data based on the quantile anal- 
ysis are consistent with those derived from the Chandra data. 

4. OPTICAL COUNTERPART 

On 2007 May 8, we observed the three Window fields 
with the Inamori Magellan Areal Camera and Spectrograph 
(IMACS) on the 6.5 m Magellan (Baade) telescope at Las 
Campanas, Chile. Under good conditions (seeing FWHM 
~ 0.5", clear sky) we obtained a dithered set of 5 point- 
ings in the f/4 configuration (15' field, 0.2"/pixel) to cover an 
18' x 18' region of BW. This provided a total exposure time 
of 500 s in each of Bessell-B, V, R, & CTIO-7 filters over the 
Chandra field. 

We processed the images using standard IRAF tasks, and 
calibrated the astrometry using the 2MASS catalogue as a ref- 
erence. The astrometric residuals on each CCD frame were 



~ 0.2". We reprojected and stacked the images using the 
SWARP 6 utility. All frames were normalized to ADU/second 
units and combined using weight-maps constructed from flat- 
fields and bad pixel masks. The initial source search and pho- 
tom etry were performed on the stacked images using SExtrac- 
tor dBertin & Arnoutsll 19961) . 

After boresight ing the initial IM ACS source list to the 
Chandra sources dZhao et alj |2005). and applying offsets of 
ARA = 0.088(51)", ADEC = 0.498(50)" to the X-ray co- 
ordinates, we found a potential counterpart in the B and V 
band images within the search area of the Chandra source. In 
Fig. [5^, the source search area is marked by the (red) circle, 
the radius of which is the 2a quadratic sum of the boresight 
error and the posi tional errors of the X-ray and optical sources 
dZhao et alj|2005l) . The potential counterpart (indicated with 
green tick marks in the figure) is not completely resolved with 
a nearby star (orange). Therefore, we generate a PSF model 
using a few bright isolated sources within ~ 30", and perform 
PSF fitting to the two sources in the B band image (where the 
candidate counterpart is brightest and its neighbor is faintest) 
using the DAOPHOT package ALLSTAR in IRAF in order to 
derive the position and magnitude of the source. The posi- 
tion of the optical source is (RA, DEC) J2000 = (18:03:54.45, 
-30:00:06.3). For the V, I and R band images, we performed 
PSF fitting with the fixed positions given by the results of the 
PSF fitting in the B band image. The instrum ental B, V and 
/ magnitudes are calibrated against OGLE-II ( Udalskij |2002l) 
stars in the field by first order polynomial fits. 

Fig-E> shows a calibrated color-magnitude diagram (CMD) 
for stars within ~ 30" of the potential counterpart (green dia- 
mond) and the neighbor (orange square). The CMD shows the 
potential counterpart is bluer than the majority of field stars. 
The calibrated apparent magnitude of the optical counterpart 
is B = 22.2(1), V = 21.7(1) and / = 20.6(1) where the errors 
are the output of the DAOPHOT. Taking the interstellar extinc- 
tion to BW to be Ay = 0-1.4, E(B-V) = 0-0.45 leads to a 
de-reddened color index (B-V)o = 0.05-0.5 for the optical 
counterpart. This value is in the typica l range for cataclysmic 
variables, see forexample|Anii(Il973), table 17.9. The num- 
ber of blue sources (24 sources with B—V < 0.7) falling in the 
error circle by chance is about 0.03. 

We estimate the ratio of the unabsorbed X-ray to optical 
flux, log(F x /F v ) = 0.63-1.12 for Ay = 0.0-1.4, where F x is in 
the 2-8 keV range. The blue color of the optical source and 
the relatively high value of log{Fx / Fy) (typically < -2.5 for 
normal stars) are consistent with the accretion nature of the 
X-ray emission from the source, indicating that the optical 
source is a good candidate for being the counterpart. 

5. DISCUSSION 

The lack of a bright optical counterpart (V < 14) rules out 
a high mass X-ray binary (HMXB) and hence a neutron star 
pulsar for the system. The relatively hard spectrum rules out a 
quiescent low mass X-ray binary or a coronal nature of the X- 
ray emission. In fact, the X-ray properties of the source show 
classic signatures of IPs such as the intrinsically hard X-ray 
spectrum with, although marginal, a possible iron emission 
line and the X-ray pulsation correlated with the internal ab- 
sorption at a period of 1028 s. The estimated luminosity of the 
system (see below) is also consistent with IPs. A likely sce- 
nario for modulating internal absorption at the spin period is 
that in a bright phase we have a relatively uninterrupted view 
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FIG. 4. — The quantile diagram of CXOPS J 180354.3-300005 as a function of modulation phase, color-coded to match the phase in Fig.[2p (see the electronic 
version). Both panels show the same data and the 10 binned data points are shown in (a). The grid in (a) is for simple power law models with T = 0, 1, 2 & 3 
and Nn22 = 0.01, 0.1, 0.4 , 1, 4 & 10 and in (b) for thermal Bremsstrahlung model with kT = 1, 2, 4, 10 keV and the same /Vh22- The arrows in (b) show the 
temporal orientation of the spectral change. There is a strong correlation between the quantiles and the modulation phase. The energy quantile E p corresponds to 
the energy below which p% of the counts are detected in the 0.3-8.0 keV band. The top-axis of the quantile diagram is labeled by the median energy (£50). 
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FIG. 5. — The finding chart in the B, V,I and R band images (left) and the color-magnitude diagram (CMD) in B vs. B-V (right). In the finding chart the (red) 
circle centered at the X-ray source (radius: 2cr of the composite errors, 0.66") indicates the search region for counterparts. The candidate counterpart and the 
neighbor star are marked by green and orange tick marks respectively (white tick marks for the print edition). The CMD compares the two marked sources with 
the field stars within 30". The potential counterpart (green diamond) is mildly bluer than the most field stars. 



of the emission region on the surface of the WD whereas in a 
faint phase we are looking at the emission region through the 
channelled accretion stream, which attenuates the X-rays. 

In the following, we estimate the likely distance of the 
source to see if the system could belong to the GB popula- 
tion, the majority of which are believed to be IPs. If we as- 
sume a typic al range of abso lute V magnitude for a CV (My 
~ 5.5-10.5) dPattersonl 19981) and the known redd ening in the 
region (A v = 1 .4 at > 3 kpc) dMarshall et al. 2006), the source 
distance should be in the range of 2-10 kpc for V ~ 22, sug- 
gesting that the source is non-local. Note if this optical source 
is not the counterpart of the X-ray source (V > 22), the dis- 
tance estimate becomes the lower bound, which reinforces the 
argument that the X-ray source is non-local. The correspond- 
ing average X-ray luminosity would be 6 x 10 31 — 10 33 ergs s" 1 
in the 2-8 keV range. According to the populati on synthesis 
model for the magnetic CVs in the GC region bv lRuiter et al.l 
(2006): iHeinke et all (120081) . this luminosity range is accept- 
able for magnetic CVs. 



We can take this a step further. Assuming the space density 
of IPs follows the stellar density, in BW we expect IPs are 
more likely found in the GB at a distance of 6-1 kpc. For 
instanc e, if we use the stellar model (model A) in lHong et al.l 
(2009), shown as (red) squares in Fig. [6] the probability that a 
given IP in the field is at 6-10 kpc is about 80% when there is 
no luminosity con straint. We can co nvert the luminosity dis- 
tribution of IPs in iRuiter et all (120061) (their standard model) 
to a probable distance distribution for the given flux of the 
source (1 .2 x 10~ 13 ergs cm -2 s" 1 in 2-8 keV), which is shown 
as (blue) diamonds in Fig. [6] Now if we combine these two 
distance-dependent constraints, we can estimate the proba- 
bility distribution of the source distance for the given X-ray 
flux (black circles). The result indicates there is about a 70% 
chance that the source is at a distance of 6-10 kpc among 
the possible range of 2-10 kpc set by the V magnitude of 
the potential counterpart (64% chance without any constraint 
from the potential counterpart). This estimate mildly favors 
the source being a bright IP in the GB X-ray population rather 
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FIG. 6. — The combined probability distribution of the source dis- 
tance (black circl es) from the stellar distribution (red squares, model A in 
Hong et al. (2009)) and the distance distrib ution based on the X -ray luminos- 
ity (blue diamonds, the standard model in Ruiter et al. 12006)) for the given 
flux of the source (1.2 X 10~ 13 ergs cirT 2 s~'). The unshaded region repre- 
sents the acceptable distance for the source based on the V magnitude of the 
potential counterpart. Each distribution is rebinned in 0.5 kpc intervals. 



than an IP outside of the GB. Note there are large uncertainties 
in the assumptions of this estimate. For instance, the statistics 
in the model of the X-ray luminosity distribution at > 10 32 



ergs s" 1 is relatively poor (see Fig. 1 in Ruit er et all (12006)). 

If indeed there are numerous IPs among the hard X-ray 
sources found in this region, long term X-ray monitoring can 
lead to a direct identification of many sources through de- 
tection of pulsations. For instance, for a 1 Ms exposure of 
the BW, we expect to detect > 500 counts from IPs with 
10 32 ergs s 1 in the GB at 8 kpc, which allows for a direct 
detection of modulations if present. A crude scaling using 
7V(> S) ~ S~ 15 predicts we should be able to identify ~ 30 
sources or more in BW from such an exposure, and at 10 32 
ergs s" 1 or below, we start to see the major part of t he luminos- 
ity dis tr ibution of the magnetic CVs according to lRuiter et al.l 
(2006); iHeinke et all (l2008l) . In addition, low extinction re- 
gions such as the Window fields are better suited for this ap- 
proach than high extinction fields such as Sgr A* since the 
dust in the high extinction fields would attenuate the appar- 
ent X-ray modulation or spectral changes dramatically, which 
would make the pulsation very hard to detect. For example, 
the spectral variation seen in Fig.|4]would not have been easily 
identifiable with the substantial external extinction (N mi ~ 6, 
a nominal value for the GC (Ba ganoff et al.l 12003b ) if the 
source was located in the GC in the Sgr A* field. 

We thank A. Ruiter for providing the X-ray luminos- 
ity distribution of IPs. This work is supported in part by 
NASA/Chandra grants GO6-7088X, GO7-8090X and G08- 
9093X. 



REFERENCES 



Allen, K. W., 1977, Astrophysical quantities., 3rd Ed., Moskva: Mir. 
Baganoff, F. K. et al. 2003, ApJ, 591, 891. 
Bertin, E. & Arnouts, S., 1996, ApJ, 117, 393. 

Cleveland, W. S., 1994, The Elements of Graphing Data, Hobart Press, 2nd 
edition. 

Drimmel, R., Cabrera-Lavers, A., & Lopez-Corredoira, M., 2003, ApJ, 409, 
205. 

Grindlay, J.E. et al. 2005, ApJ, 635, 907 
Grindlay, J.E. et al. 2009, in preparation. 
Freeman, P.E. et al. 2002, ApJS, 138, 185. 
Gehrels, N. et al. 1986, ApJ, 303, 336. 
Heinke, C. O., et al., 2008 AIPC, 1010, 136. 

Hong, J, Schlegel, E. M. & Grindlay, J. E., 2004, ApJ, 614, 508. (H04) 

Hong, J, et al. 2005, ApJ, 635, 907. (H05) 

Hong, J, et al. 2009, submitted to ApJ. 

Hong, J, et al. 2009b, in preparation. 

Laycock, S. et al. 2005, ApJL, 634, 53 (L05) 

Marshall, D. J. et al. 2006, A&A, 453, 635 

Muno, M. P. et al. 2003, ApJ, 589, 225 (M03) 

Muno, M. P. et al. 2003a, ApJ, 589, 225 



Muno, M. P. et al. 2003b, ApJ, 599, 465 
Muno, M. P. et al. 2004, ApJ, 613, 1 179 (M04) 

Norton, A. J., Wynn, G. A. & Somerscales, R. V., 2004 ApJ, 614, 349 

Patterson, J., 1998, PASJ, 110, 1132. 

Predehl, P. & Schmitt, J. H. M. M., 1995, A&A, 293, 889. 

Ramsay, M. & Cropper, M. 2004, MNRAS, 347, 497. 

Ritter, H. & Kolb, U., 2003, ApJ, 404, 301 

Ruiter, A., Belczynski, K. & Harrison, T. 2006, ApJL, 640, 167. 

Scargle, J. D., 1982, ApJ, 263, 835. 

Schlegel, D., Finkbeiner, D. & Davis, M., 1998, ApJ, 500, 525. 

Sumi, T., 2004, MNRAS, 349, 193. 

Udalski, T, 2002, ACTA Astronomica, 52, 217. 

Wang, Q. D. , Gotthelf, E. V. & Lang, C. C, 2002, Nature, 415, 148. 

Watson, M. G, et al., 2009, A&A, 493, 339. 

van den Berg, M. et al. 2006, ApJL, 135, 647 

van den Berg, M. et al. 2009, submitted to ApJ. 

Zhao, P. et al. 2005, ApJS, 161, 429. 



